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ABSTRACT 

We investigate the properties of quiescent and star-forming galaxy populations to z ~ 2 with purely photo- 
metric data, employing a novel rest-frame color selection technique. From the UKIDSS Ultra-Deep Survey 
Data Release 1, with matched optical and mid-IR photometry taken from the Subaru XMM Deep Survey and 
Spitzer Wide-Area Infrared Extragalactic Survey respectively, we construct a Zf-selected galaxy catalog and 
calculate photometric redshifts. Excluding stars, objects with uncertain z p h t solutions, those that fall in bad or 
incomplete survey regions, and those for which reliable rest-frame colors could not be derived, 30108 galax- 
ies with K < 22.4 (AB) and z p h t < 2.5 remain. The galaxies in this sample are found to occupy two distinct 
populations in the rest-frame U — V vs. V — J color space: a clump of red, quiescent galaxies (analogous to 
the red sequence) and a track of star-forming galaxies extending from blue to red U — V colors. This bimodal 
behavior is seen up to z ~ 2. Due to a combination of measurement errors and passive evolution, the color- 
color diagram is not suitable to distinguish the galaxy bimodality at z > 2 for this sample, but we show that 
MIPS 24/im data suggest that a significant population of quiescent galaxies exists even at these higher redshifts. 
At z = 1-2, the most luminous objects in the sample are divided roughly equally between star-forming and 
quiescent galaxies, while at lower redshifts most of the brightest galaxies are quiescent. Moreover, quiescent 
galaxies at these redshifts are clustered more strongly than those actively forming stars, indicating that galaxies 
with early-quenched star formation may occupy more massive host dark matter halos. This suggests that the 
end of star formation is associated with, and perhaps brought about by, a mechanism related to halo mass. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — infrared: 
galaxies 
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1. INTRODUCTION 

The stellar mass in the local universe lies largely within 
two distinct classes of galaxies: actively star-forming spi- 
rals, and massive ellipticals with evolved stellar popula- 
tions and little current star formation. These populations 
show strongly bimodal behavior in a n umber of measure d 
and derived quantities, including color (Ba ldry et al.ll2004h . 
the 4000A-break s trength £> n (4000) (an i ndicator of stel- 
lar popu lation age; Kauffman n et al.l 120031) . and clustering 
strength dBudavari et al.ll2.Q03b . Relatively low-mass galax- 
ies with high star-formation rates have also been found at 
high redshifts (e.g ., z > 3 Lyman-break galaxies or LBGs; 
ISteidel et al.lll996h . Surprisingly, the population of massive, 
"dea d" galaxies appears to persist to high redshift as well 
(e.g. lLabbe et al.ll200H IPaddi et al.ll200l iKriek et al.ll2006l 
2008b), even though the universe was only a few Gyr old at 
that point. Since the stellar populations of evolved galaxies 
often have ages that are a significant fraction of the age of the 
universe at these redshifts, the implication is that they must 
have formed (and their star formation ceased) far earlier than 
anticipated by standard hierarchical structure formation mod- 
els. 

Investigating the evolution of star formation activity and the 
assembly of stellar mass over the past ~ 12 Gyr are there- 
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fore two key avenues to understanding how the present galaxy 
population came to be, and large samples of galaxies during 
this peak formation epoch (z ~ 2) are needed for this. As spec- 
troscopy of faint objects is extremely expensive with regard 
to telescope time, several highly efficient broadband selec- 
tion techniques have been designed to weed out low-redshift 
interlopers. In addition to the aforementioned LBG selec- 
tion, which primarily finds unobscured star-forming galaxies 
at z > 3, the distant re d galaxy (DRG) c riterion (J—K > 2.3 
(Vega) or > 1 .34 (AB); lFranx et alj2003l) i nstead tends to find 
the m ost massive galaxies at 2 < z < 3 (Ivan Dokkum et all 
2006), about h alf of which show si gns of heavily-obscured 
star formation dPapovich et alj|2006l). O ther methods, such as 
the BzK selection of Dad di et all (120041) . are adept at selecting 
z ~ 2 galaxies over a wide range of masses and star-formation 
rates. 

Each of these techniques has its advantages and disadvan- 
tages, but one point has become clear: to effectively study the 
mass evolution of high-redshift galaxies, deep near-infrared 
data are crucial. The observed JHK bands trace rest-frame 
optical light at z ~ 2 - 3 and provide a more reliable indica- 
tor of stellar mass than the rest-frame UV (observed optical); 
indeed, DRGs themselves appear to represent an important 
fraction, if not the majority, of stellar ma ss at high redshifts 
dRudnick et aUl2006t iGrazian et alJl2007l) . Such galaxies are 
typically faint in the observed optical bands, either from old 
stellar populations or dust obscuration, and are often misse d 
by purely optical selection techniques (Qua dri et alJ l2007b). 
Selecting galaxies from deep near-infrared fields is therefore 
likely to yield a more complete picture of the stellar mass at 
high redshift than similar studies at optical wavelengths. 

Until recently, the small sizes (at most a few square arcmin- 
utes) and relatively low efficiencies of the available infrared 
detectors made large, deep near-IR surveys impractical. Such 
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projects therefore typically followed the "pencil-beam" ap- 
proach, surveying a single frame to high sens itivity (e.g., 
the F aint Infrared Extragalactic Survey, FIRES; iFranx et alj 
2003), or alternatively observing som ewhat larger areas at th e 
expense of depth (e.g., GOODS-S: iGiavalisco et al.ll2004l) . 
While such surveys have revealed a great deal about the 
high-redshift universe, the necessarily small survey volumes 
proved problematic for statistical studies of galaxy popula- 
tions (as a result of both cosmic variance and small-number 
statistics). 

Large-format near-IR detectors on 4m-class telescopes 
(such as ISPI, WFCAM, and the upcoming VISTA cam- 
era) now provide the ability to conduct deep surveys over 
wide areas, and a number of past, present, and planned 
projects take advantage of this capability - examples in- 
clude the Multiwave l ength Survey by Yal e-Chile (MUSYC; 
iGawiser et al.1 l2006t lOuadri et all l2007bh. the UKIRT In - 
frared Deep Sky Survey (UKIDSsT lLawrence et ail 120071) . 
and ULTRA-VISTA. Until this latter survey is complete, 
the UKIDSS Ultra-Deep Survey (UDS 5 ) is the premier 
near-IR dataset in terms of depth and area. Overlapping 
deep optica l imaging from the Subaru-XMM Deep Sur- 
vey (SXDS; ISekiguchi et al.ll2004l) and shallow mid-IR data 
from the Spitzer Wid e-Area Infrared Extragalactic (SWIRE; 
lLonsdale et all 120031) survey provide complementary data 
over a broad range of wavelengths. 

Here we derive a multiband ^-selected galaxy catalog from 
the UDS, SXDS, and SWIRE data, calculating photometric 
redshifts and rest-frame colors for all objects in the overlap- 
ping survey area. Rather than rely on a specific color selection 
technique, we instead define samples of galaxies at various 
redshifts directly from the computed z p h t values, and further 
subdivide the samples into star-forming and quiescent galax- 
ies based on their rest-frame colors. The color evolution with 
redshift, and the clustering of 1 < z p h t < 2 galaxies, are then 
determined. Readers who are primarily interested in these sci- 
ence results can find this discussion beginning in $6] 

We present an overview of the survey data in §|2] and the 
process of preparing and matching the datasets to each other 
in 5(3] The extraction of the catalog is then discussed in 
and the derivation of photometric redshifts and rest-frame 
colors in $5] $6] describes the bimodality of galaxies in rest- 
frame color space and the separation into quiescent and star- 
forming samples, and the clustering properties of these two 
populations are discussed in $7] Further discussion, caveats, 
and a summary of the results can then be found in $8] and 
jfTOl respectivelv. 

A cosmology with Hq = 70km s -1 Mpc -1 , fi m = 0.3, and 
JIa = 0.7 is assumed throughout, and magnitudes are quoted 
in the AB system unless otherwise noted. 

2. DATA CHARACTERISTICS 

In this analysis we employ the reduced J- and Zf-band 
UDS mosaics p rovided as a subset of the UKIDSS Data Re- 
lease 1 (DR1; IWarren et al.ll2007al) ; fl-band data were not 
available in thi s relea se. The UKIDSS project, defined in 
lLawrence et al.1 (j2007l). uses the UKIRT Wide Field Camera 
(WFCAM ICasali et all |2007h a nd a photometric system de- 
scribed in iHewett et al.l d2006i) ; the pipeline processing and 
scie nce archive are desc ribed in Irwin et al. (in preparation) 
and lHambly et al.l(l2008l) respectively. The UDS field consists 
of a single repeatedly-observed UKIDSS survey tile (in turn 

s |http : / /www . nottingham . ac ■ uk/astronomy/UDS/| 




FIG. 1 .— UKIDSS K-band image of the UDS field, with SXDS (solid line) 
and SWIRE (dashed line) coverage overplotted. The image is about 52' on a 
side; north is up and east is to the left. 

comprising four subsequent WFCAM observations arranged 
in a square to produce uniform coverage), with a total area 
of 0.77 deg 2 and reaching nominal depths of K < 23.45 and 
J < 23.55 (AB; 5a point-source threshold in a 2" aperture). 
The 0."4 pixel scale of WFCAM does not sufficiently sample 
the A"-band PSF of ~ 0."7, so the individual observations are 
offset by fractions of a pixel and the combined mosaics regrid- 
ded by a factor of 3 to a pixel scale of 0"134. No te that while 
Data R elease 2 has since been made available (IWarren et alJ 
2007b), this release contains identical UDS mosaics to DR1; 
only data in the shallower UKIDSS surveys were updated. 

Optical imaging of most of the UDS field is provided by 
the Subaru-XMM Deep Survey (SXDS) beta release 6 , a deep 
wide-field survey employing the Subaru telesc ope and Subaru 
Prime Focus Camera (Suprime-Cam or SUP; iMivazaki et alj 
120021) . This imager contains ten 2048 x 4096 pixel CCDs 
arranged in a 5 x 2 mosaic, providing a rectangular field- 
of-view of roughly 34' x 27'. The SXDS then consists of 
five overlapping SUP fields arranged in a cross-shaped pat- 
tern (three adjacent fields in the north-south direction, and 
two on the eastern and western flanks with the camera rotated 
90°; see Figure[T|i. Reported 5a depths (in a 2" aperture) are 
B < 27.5, R < 26.8, if < 26.5, and z' < 25 A, though the sen- 
sitivity varies by ~ 0.2 between the five fields (with a ~ 0.6 
mag variation in z); our independent sensitivity estimates are 
comparable to the reported values. 

The SWIRE survey covers a total of — 49 deg 2 with IRAC 
and MIPS in six subfields, one of which (the XMM-LSS 
field) includes nearly the entire UDS/SXDS area. The orien- 
tations of the IRAC/MIPS fields are at an angle to the UDS; 
therefore, only a small triangular region in the southwest cor- 
ner of the UDS has no SWIRE coverage (shown in FigureQ]). 
SWIRE reaches nominal 5a point-source depths of 3.7, 5.4, 
48, and 37.8 ^Jy in the 3.6, 4.5, 5.8, and 8.0 ^m IRAC bands 
respectively, and 230 /iJy in MIPS 24 ^m. The IRAC mo- 
saics of this field are divided into 16 separate subimages, with 
four of these having some overlap with the UDS field; the 
MIPS 24 /im image of the field is available as a single file. 
The relevant subimages of all four IRAC bands, as well as 

6 http : //www . nao j . org/ Science /SubaruPro ject /SXDS/ 
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the full 24 fim MIPS image, were obtained from the pub- 
lic SWIRE data website 7 . Due to their far broader point- 
response functions and lower sensitivities, we do not analyze 
the MIPS 70 /im or 160 /im data here; similarly, the 5.8/im 
and 8/im data are excluded due to their order-of-magnitude 
lower sensitivity compared to the first two IRAC bands. Fur- 
thermore, since most objects seen in the UDS K-band are not 
individually detected in the SWIRE 24/im image, the 24/im 
data are therefore not included in our final ^-selected cat- 
alog (but a stacking analysis of these data are discussed in 
36.4b . A summary of the data employed in the production of 
our /f-selected catalog, inc ludin g areas and measured limit- 
ing depths (as discussed in 34.51 ) is given in Table[T] 

3. DATA PREPARATION 
3.1. Image registration 

The simplest way to achieve consistent flux measurements 
across multiple bands with arbitrary PSFs is to transform the 
images such that they are all on the same coordinate sys- 
tem and pixel scale. The pixel scale of the UDS images is 
0."134; however, all of these images have PSF FWHMs larger 
than 0."7 (the UDS K-band seeing). To reduce file sizes 
and processing time, we thus adopt a standard pixel scale of 
0."202, the native scale of the SXDS images. This is some- 
what coarser than the original regridded UDS pixel scale, but 
still sufficiently samples all PSFs. 

The SWarp software package 8 (B ertin et alj2 002) was used 
to crop, resample, and transform these sets of images to a 
common coordinate system before combining the subimages 
(in the cases of the IRAC and SXDS data) into single mo- 
saics. Since the UDS K image has the best seeing and is 
the band from which objects are detected, this was taken as 
the "standard" field to which all other images were regis- 
tered. The J and K data were resampled to the SXDS pixel 
scale and cropped to remove noisy and missing data regions 
near the edges, resulting in 52'A by 53 ^0 images. Likewise, 
the SWIRE MIPS 24 fim image and the IRAC subimages 
that cover the UDS field were transformed, resampled, and 
cropped to cover the same area. However, since the five SXDS 
pointings taken with each filter can exhibit different seeing 
and geometric distortions, those images were PSF-matched 
(see 33.2b and astrometrically corrected prior to combining 
with SWarp. 

As SWarp transformed the images based solely on the as- 
trometric information in the image FITS headers, some dis- 
tortions are likely to remain (particularly between data from 
differing instruments). The measured positions of stars were 
thus used to correct th ese distortions. Sourc e Extractor ver- 
sion 2.5.0 (SExtracton lBertin & Arnoutsll 19961) was first used 
to find the positions of bright objects in the K image and mea- 
sure their approximate J and K fluxes. 1179 stars were se- 
lected on the basis of their blue IR colors (/— K < 0) and 
bright but unsaturated K-band fluxes (17.4 < K < 21 .4). The 
radial profiles of these objects were fitted with Gaussians to 
determine their centers and PSF widths in the K band; all but 
a negligible fraction (< l%)had FWHMs comparable to the 
nominal image resolution, indicating that this simple color cut 
indeed primarily selects point sources. Next, a centroid algo- 
rithm was employed to determine the positions of these stars 

7 |http: // swire ■ ipac ■ caltech . edu/swire/astronomers/] 
|da ta_access.html| 

° See |httpT77terapix . iap ■ f r/rubrique ■ php?id_rubrique=4 
for further information 



in each of the BRi'z'JK/3.6/4.5 images. Finally, the pixel 
positions of these stars along with the IRAF tasks geomap 
and geotran were used to transform the SXDS, UDS /, 
and IRAC images to the same coordinate system as the UDS 
K mosaic. To check the agreement between the images, the 
pixel coordinates of these stars were again measured in the 
transformed images and found to agree well, e.g. in the B mo- 
saic the image coordinates agree with those in the K image to 
within 0.5 pixel (0."1) in 75% of cases, with a median abso- 
lute deviation of 0."06. Similar offsets are seen in the trans- 
formed IRAC 3.6 and 4.5/im images (in addition, automatic 
astrometric corre ction s are performed in the IRAC deblend- 
ing described in 34.21 ). Such astrometric deviations typically 
translate to systematic flux errors of at most ~ 1%. 

3.2. IR/Optical PSF matching 

The observed seeing width varies significantly between the 
different bands: point sources in the UDS have FWHMs of 
approximately 0"7 in K and 0."8 in J, while in SXDS they 
range between 0."7 — 0."9 depending on the field and band. 
Additional higher-order variations are likely to be present in 
the PSFs due to differences in the two instruments' optics. Di- 
rectly performing aperture photometry on these images would 
lead to flux offsets in the different bands (with less flux falling 
inside the aperture for the poorer-seeing images), and thus 
systematic color errors. 

To mitigate this, all images were smoothed to the same PSF 
before fluxes were measured. Although this can be accom- 
plished by convolving each with a simple Gaussian or higher- 
order kernel, the multitude of stars within these images can 
themselves be used to generate empirical kernels that con- 
tain all structural information about the PSFs. This is accom- 
plished by deconvolving a low-resolution ("reference") PSF 
with one constructed from a higher-resolution image: 

K = P l - ef (g>- 1 P hiKS (1) 

where P re f and Phires are the low-resolution and high- 
resolution PSFs respectively. The resulting kernel can then 
be convolved with the high-resolution image in its entirety to 
bring it to the same PSF as the reference image. 

ImMch(x,y) = /hires (x,y)®K (2) 

Since the deconvolution step strongly magnifies any intrinsic 
noise in the PSFs, a large number of bright (but unsaturated) 
stars are necessary in each band to produce PSFs with suffi- 
ciently low noise levels. 

For the J and K images, the PSFs were constructed from 
the same stars that were selected for astrometric matching, 
keeping the 200 brightest in each band. However, many of 
these stars are saturated in the SXDS images due to the much 
longer exposure times employed in this survey. As the effects 
of saturation are far more detrimental for the PSF shape than 
in the astrometric correction, a different method for select- 
ing stars in the SXDS images is necessary. For these images, 
the IRAF psf measure task was used to estimate widths of 
unsaturated objects in the image based on a Gaussian fit. In 
a histogram of the widths of these objects, there is an obvi- 
ous cutoff on the low-FWHM end accompanied by a sharp 
peak and an overlying broad distribution extending to larger 
FWHMs. Objects in the narrow peak were presumed to be 
point sources to within the image resolution, and the brightest 
of these were selected as stars in each optical image. 
9 Square postage stamps (21 pixels or 4 ."2 on a side) of the 
stars selected in each image were then created, and these in 
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FIG. 2. — Average growth curves of the six optical/IR bands divided by the 
growth curve of the worst image (SXDS-West z-band, with a PSF FWHM 
~ 0"9). Relative growth curves before (dashed) and after (solid with points) 
convolution to the broadest PSF are shown. The vertical dotted line shows 
the "color" aperture size of 1 ''75. Typical systematic flux errors due to PSF 
matching at this aperture size are < 2%. 

turn were median-combined to create 22 empirical PSFs (one 
for each of the 20 SXDS subimages, plus one for each UDS 
mosaic). Enclosed flux as a function of aperture size (growth 
curves) were then created for these PSFs. The "slowest- 
growing" growth curve, indicating the broadest PSF, is that 
of the western SXDS z-band image with a FWHM of ~ 0"9. 
Thus, we take this as the reference PSF and match the other 19 
SXDS fields and two UDS images to it. The average growth 
curves for point sources in the smoothed BRi'z'JK images, ex- 
pressed as a fraction of the SXDS-West z' (reference) growth 
curve, before and after smoothing, are shown in Figure [2] At 
the aperture sizes used for flux measurements in the follow- 
ing section (> 1 ."75 or ~ 9 pixels in diameter), the PSFs are 
matched to within 2%. 

The point-response functions of the IRAC 3.6 /im and 
4.5 /.tm images are roughly 1.7x as broad as in the worst 
SXDS image (~ 1"6 for both); they also contain significant 
non-Gaussian structure resulting from point source diffrac- 
tion. While it would have been simpler in principle to smooth 
the BRi'z'JK images to the IRAC PSF shape, this would 
substantially reduce the detection efficiency in K while (for 
most high-redshift sources) providing only weak detections 
in IRAC to begin with. Additionally, blending from nearby 
sources is more likely to adversely affect simple aperture pho- 
tometry in the IRAC images. PSF-matching and photometry 
are thus performed on these data using a different technique, 
described in 34.21 

4. CATALOG EXTRACTION 

4.1. Source detection and photometry 

To measure fluxes of sources in the BRi'z'JK images we 
used SExtractor in dual-image mode, whereby one image is 
used to detect sources and aperture photometry is performed 
on another image registered to the identical coordinate sys- 
tem. The unsmoothed (pre-PSF matching) K mosaic was 
used as the detection image, and fluxes were measured in each 
of the six PSF-matched IR/optical mosaics. The sensitivity 
across the K image is not perfectly uniform due to overlap- 
ping exposures and array efficiency variations; such nonuni- 



FlG. 3. — Top panel: Median growth curves of point sources in the PSF- 
matched BRi'z'JK images; Bottom panel: Signal-to-noise in the six bands, 
defined as the growth curve divided by the empty aperture noise function in 
each band, normalized to the peak SN in each band. The solid vertical line 
denotes the color aperture size of 1 "75. 

formity can lead to corresponding variations in the number of 
detected galaxies, and thus to spurious clustering signals. We 
used the UKIDSS-supplied confidence map to construct an 
RMS map of the K image, and this map was applied as the 
weight image within SExtractor in order to detect sources in 
an effectively noise-equalized and uniform manner. Several 
initial SExtractor runs were performed with varying detection 
parameters and the results checked by eye; the final choice of 
parameters appe ared to find all faint objects with few spurious 
detections (see 34.51 >. 

Total fluxes were determined usin g a flexibl e elliptical aper- 
ture (SExtractor's AUTO aperture ; iRronl 1980). To ensure that 
color measurements are consistent, fluxes in a fixed circular 
aperture were also measured. The optimal size of this "color" 
aperture is subject to two primary considerations: while a 
larger aperture encloses a greater fraction of a given object's 
flux, it also suffers from higher uncertainty due to background 
fluctuations. To determine the best aperture size, we divided 
the growth curve of each image by its noise function (i.e., 
photometric uncertainty versus aperture size, described in de- 
tail in 34.4b . This is effectively equivalent to the signal-to- 
noise (SN) of each image as a function of aperture size, and is 
shown in the lower panel of Figure [3] 

The /T-band SN peaks at a diameter of ~ 4 pixels or 0."8, 
which is lower than the ideal aperture size theoretically ex- 
pected with uncorrelated noise (~ 1 ."0), but consistent with 
the 1 . 1 - 1 .4 x FWHM opti mal aperture found for the MUS YC 
data dOuadri et a l. 2007b). The growth curves shown in the 
top panel of Figure [3] indicate that an aperture of this size 
misses about 70% of a point source's flux. Such a small aper- 
ture is highly susceptible to systematic errors from imperfect 
PSF matching and astrometric offsets. Additionally, larger 
apertures appear to produce slightly bluer J- K colors (on the 
order of 0.05 magnitudes between 1" and 2" apertures), per- 
haps due to intrinsic galaxy color gradients and/or imperfect 
image matching. To ensure more accurate colors of non-point 
sources, as a compromise we adopt a somewhat larger "color" 
aperture of 1 ."75 diameter. This more than doubles the en- 
closed flux while only decreasing the SN by ~ 30% from the 
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FIG. 4. — Distribution of total K magnitudes for objects in the UDS K- 
selected catalog detected at > 5<r significance in the "color" (1"75 diameter) 
aperture black histogram. The vertical black line indicates the mean 5a detec- 
tion limit of K = 23.57. O verplotted are raw number co unts from the FIRES 
MS 1054-03 cluster field I Forster Schreiber et al. 2006, red histogram) and 
the CDFS-GOODS jf-selected catalog jWuvts et aU2003 , green histogram). 



optimal value. 

In cases where the AUTO aperture is smaller than 1 ''75 
(primarily occurring for faint, point-like sources), the color 
aperture flux is substituted for the total flux. All total fluxes 
are corrected for flux falling outside the aperture assuming a 
minimal (point-like) source flux distribution. The SExtractor 
photometry was verified by comparing the z' total flux values 
with those of objects in spatially coincident HST/ACS im- 
ages. 9 The only useful ACS data in the archive were taken 
with the F850LP filter, which is approximately equivalent to 
SXDS z'; fluxes of objects in these ACS data were consistent 
with those derived from SXDS. The photometric zeropoints 
of the SXDS z' data thus appear to be properly calibrated. 
Since the same calibration techniques and standard star fields 
were used by the SXDS team to determine the BRi' image 
zeropoints, we assume the calibrations in these bands are ac- 
curate as well (though the lack of available verification data 
should be kept in mind). Figure [4] shows the distribution of K 
magnitudes for all objects detected by SExtractor, with num- 
ber counts overplot ted from two other A"-select ed samples: 
FIRES M S 1054-03 dForster Schreiber et aU2006h and CDFS- 
GOODS dWuvts et alj|2008l) . The number counts of the three 
surveys are similar to the UDS limiting magnitude, though the 
enhanced number of faint objects in the FIRES field (due to 
the presence of the z = 0.83 cluster MS 1054-03) is evident. 

4.2. IRAC PSF-matching and photometry 

As previously mentioned, the optical and NIR images were 
smoothed to the seeing of the SXDS-West z'-band data rather 
than the much broader point-response function of the IRAC 
data. The direct aperture photometry from SExtractor would 
therefore not measure IRAC fluxes that are properly matched 
to the other bands. Moreover, the lower-resolution IRAC im- 
ages are more prone to blending than the other bands, which 
can introduce significant systematic err ors in the measured 
colors. We thus employed the method of Labbe et al.l d2006l) 

9 Found using the MAST archive, [http : / / archive . stsci . edu| 



to correct this, summarized below. 

Sources that are bright in K are also typically bright in 
the 3.6/4.5 /im bands; thus, the K image is used as a high- 
resolution template to deblend the IRAC images. Convolution 
kernels used to transform from the K to the IRAC PSFs are 
constructed from bright and unsaturated sources (computed 
by fitting a series of Gaussian-weighted Hermite functions to 
the Fourier transforms of the objects in IRAC and K), and a 
smoothed map of the kernel coefficients as a function of im- 
age position is then produced. For a given object detected 
in K, neighboring sources are fitted and subtracted out using 
the local convolution kernel (derived from the aforementioned 
smoothed coefficient map) and the SExtractor-derived seg- 
mentation map. This results in a "cleaned" IRAC subimage 
containing only the source in question; standard photometry 
with a 3" diameter aperture 10 is then performed on the source, 
and this process is repeated for each object in the /^-selected 
catalog. Through a visual inspection of the IRAC residual im- 
age with all sources subtracted, we conclude that this method 
effectively removes contaminating so urces (for an illust ration 
of this technique, see also Figure 1 in lWuvts et alJl2007D . 

Since the IRAC photometry from this procedure is per- 
formed using a larger aperture on images with significantly 
broader PSFs than the matched BRi'z'JK images, these fluxes 
must be corrected in order to obtain accurate IRAC-IR/optical 
colors. This is accomplished by multiplying the IRAC flux 
measured in a 3" aperture by a correction factor to obtain a 
"matched" IRAC flux: 

/k(1"75) 



/iRAC. match = /lRAc(3 ) X 



(3) 



.conv (3") 

where /k.conv(3") is the flux of the object measured in the K- 
band image after smoothing it to the IRAC resolution, and 
/k is the K-band flux measured in the fixed "color" aperture 
(see 94.lt . This correction implicitly assumes that the ratio 
of fluxes between the two aperture sizes is the same for K 
and IRAC, i.e. color gradients between 1 ."75 and 3" are in- 
significant. While this technique could introduce systematic 
errors in measurements of extended sources with strong color 
gradients, at the redshifts primarily considered here (z > 1) 
most galaxies have small angular sizes and are not likely to 
be adversely affected. 

4.3. Bad pixels 

All three sets of images (UDS, SXDS, SWIRE) contained 
regions where either no data were available or the derived 
photometry was unreliable. Thus, even though SExtractor 
finds objects in these regions, their fluxes should be consid- 
ered unreliable. To account for this we created bad-pixel 
maps for each image set; if an object's "color" aperture con- 
tains any such bad pixels, we set a flag indicating as such in 
the catalog. The bad pixel maps created for each survey are 
mostly a result of the following effects: 

1. UDS JK:False sources and large negative residuals. In 
rows and columns containing bright stars, there are of- 
ten significant positive and negative images (repeated at 
regular intervals) as a result of crosstalk between WF- 
CAM detector nodes. These "false positive" images, 
while obviously not true astrophysical sources upon vi- 
sual inspection, are nonetheless picked up as real ob- 
jects by SExtractor. To make matters worse, many such 

10 As with the UDS and SXDS images, this aperture size was chosen as 
a compromise between maximal signal-to-noise and enclosing a sufficient 
fraction of a point source's flux. 
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"false positives" in K have negative apparent flux in 
J, causing these to appear as very red objects. They 
can thus mimic high-redshift near-IR galaxy colors and 
may be an important contaminant in our catalog. Fortu- 
nately, adjacent to most of the false positive images in 
K, some strongly negative residuals are also typically 
seen. We thus searched for such strong K residuals 
in the image, removing any object falling on or within 
10 pixels of any such artifact as bad. As a secondary 
check, we searched for objects which were detected in 
K but undetected in the sum of the B+R+i'+z' mosaics. 
This cut finds an additional ~ 1900 possible bad spots. 
A visual inspection confirms that about 80% of these 
are indeed artifacts (including diffraction spikes from 
bright stars and meteor streaks in K) that were unde- 
tected by the aforementioned "negative residual" tech- 
nique, and most of the remaining 20% are extremely 
faint and likely to be spurious detections. However, 
since some of these optically-nondetected objects are 
real, they have been included in the catalog with a flag 
indicating that they are likely to be bad. 

2. SXDS BRi'z 1 '- Non-covered regions and bright stars. 
As shown in Figure [2 the corners of the UDS fields are 
not covered by the SXDS imaging. Additionally, bright 
stars in SXDS are surrounded by concentric halos and 
extended readout streaks. The Subaru team has pro- 
vided files which define the regions affected by these 
bright stars; a visual inspection of the SXDS frames 
confirms that these regions accurately describe the af- 
fected areas. For each of the five pointings we created 
a map such that pixels falling within these regions had 
a value of 1, and zero otherwise. Using SWarp, these 
bad pixel masks were then transformed and combined 
into a single mosaic on the same coordinate system as 
the other images (assigning a value of 1 to the corners 
with no SXDS coverage as well). 

3. SWIRE: One non-covered region. The SWIRE survey 
spans an area many times larger than the UDS field, 
with nearly uniform coverage over the UDS/SXDS area 
except for a small triangular region in the southwest 
corner of the field. As most of this triangle is also out- 
side the SXDS coverage area, reliable IRAC flux mea- 
surements (or upper limits) are available for essentially 
all objects with SXDS photometry. Objects with coor- 
dinates falling in this missing SWIRE piece are flagged 
in the catalog. 

It should be noted that, while these techniques were effec- 
tive at automatically finding many false detections and areas 
of unreliable photometry, some undetected artifacts are likely 
to remain in the full catalog. For the analysis described herein 
we therefore visually inspected subsets of our samples to en- 
sure that such artifacts do not significantly affect our results. 

4.4. Photometric Error Determination 

The uncertainty of a flux determined through aperture pho- 
tometry contains contributions from the intrinsic Poisson er- 
ror in the source counts and from background noise. While 
the first (automatically calculated by SExtractor) is straight- 
forward to determine and is the dominant source of uncer- 
tainty for bright objects, the background uncertainty is impor- 
tant for fainter sources. For perfectly random (uncorrelated) 
background noise the additional photometric uncertainty is 



TABLE 1 

UDS/SXDS/SWIRE DATA CHARACTERISTICS 



Survey 1 


Band 2 


Area 3 


4 

>"5<T 






(deg 2 ) 


(AB) 


SXDS 


B 


0.70 


27.65 


SXDS 


R 


0.70 


27.05 


SXDS 


i' 


0.70 


26.82 


SXDS 


z> 


0.70 


25.53 


UDS 


.1 


0.77 


23.93 


UDS 


K 


0.77 


23.57 


SWIRE 


3.6/jm 


0.76 


22.25 


SWIRE 


4.5/im 


0.76 


21.53 



1 SXDS — Subaru-XMM Deep Survey; UDS— UKIDSS 
Ultra-Deep Survey; SWIRE — Spitzer Wide-Area Extragalac- 
tic Survey. 

2 The photometric systems are: Johnson-Cousins (BIR); SDSS 
(i'/z'); Mauna Kea (J/K). 

3 Survey areas are in square degrees overlapping the UDS. 

4 Limiting magnitudes (5a) are calculated from the average 
noise properties of each image using the 1 "75 color aperture, 
and scaled up to the total object flux using a minimal aperture 
correction factor (i.e., assuming a point source) derived from 
the growth curves in Figure [5] 

simply (Taper ~ pVN, where N is the number of pixels con- 
tained within the aperture and a is the pixel-to-pixel RMS. 
On the other extreme, if the noise is perfectly correlated, then 
Caper ~ crN . Typical images will fall somewhere between 
these two cases, exhibiting a general form of er aper ~ a(\/N)P. 

We estim ated the backgroun d unc ertainty follow i ng th e 
method of lLabbe et all (120031) and iGawiser et all (120061) . 
First, to estimate this power-law index, several hundred aper- 
tures of a given size were randomly placed on empty parts 
of the images (i.e., those containing no emission detected by 
SExtractor). The RMS variation cr aper was found by fitting 
a Gaussian to the resulting histogram of aperture fluxes, and 
this process was repeated for aperture diameters from 2 to 40 
pixels (0"4-8"). We then fit a general power-law function 
to the cr a per vs. \/N curve to determine the average value of 
j3 over each image; these (3 parameters were near 1.5 for all 
images (ranging from (3=1 .37 for J to (3 = 1 .58 for ;')■ 

Ideally, the power-law normalization should simply be 
equal to the local pixel-to-pixel RMS a. However, in real 
images some deviation from this is possible. To accurately de- 
termine this relation, we performed the aforementioned pro- 
cedure over ~ 300 random, small (2000 x 2000 pixel) subre- 
gions of each image, this time fitting a function of the form 
c aper = aaiy/N) 13 and fixing (3 to the average value measured 
for each mosaic. The value of a for the pixels contained 
within the apertures was also calculated, and a linear fit per- 
formed to determine the constant of proportionality a. For 
the six images a averages about 1.3 (ranging between 1.22 
and 1.40), reasonably near the theoretically-expected value 
of 1. For each object in the catalog, the local value of a 
and the above formula was used to estimate the background 
uncertainty, and this was combined in quadrature with the 
SExtractor-derived Poisson noise to compute the total flux 
uncertainty. 

4.5. Limiting magnitudes and completeness 

From the noise estimates we calculated approximate So- 
limiting magnitudes in the 1 "75 color aperture for each band, 
listed in Table Q] These are based on the ape rture noise mea- 
surements in each image as described in 34.41: the actual limit- 
ing magnitudes vary as a function of position in each mosaic, 
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FIG. 5. — Point-source completeness of the ^-selected catalog, for the 
average mosaic depth (solid line) and the least sensitive part of the mosaic 
(dotted line). 

particularly in the SXDS mosaics. The magnitudes listed in 
the table have been corrected for flux falling outside the color 
aperture (assuming a point source), as determined from the 
growth curves shown in Figure [3] Note that these limiting 
magnitudes are slightly fainter than those reported by the sur- 
vey teams; this is due to our use of a smaller color aperture 
(l."75 vs. 2"). 

Completeness was tested by placing simulated point 
sources in a 5000 x 5000 pixel subregion of the unconvolved 
UDS K image. This subregion sits exactly in the center of 
the full mosaic, and was chosen to contain roughly equal 
contributions from the four WFCAM arrays (which individ- 
ually have slightly different sensitivities), thereby providing 
a representative sample of the full mosaic. The K PSF gen- 
erated for the PSF matching step was used to create the sim- 
ulated stars. Approximately 1700 such stars, scaled to mag- 
nitudes between K= 18-23, were placed in the image on 
a semi-random grid designed to avoid other real or simu- 
lated sources. SExtractor was then run on the image with 
the point sources included (using identical detection param- 
eters to those in 34.11 ). and the number of detected fake stars 
as a function of magnitude recorded. At the mean 5a Zf-band 
detection limit of K = 23.57, the survey is 95% complete. 

This process was repeated for one of the lowest-sensitivity 
regions in the mosaic, the area of a single WFCAM chip in the 
northwest quadrant of the image. Even in this worst case the 
catalog is 100% complete for K < 23.2 and 95% complete for 
K < 23.5 (see Figure 0. It should be noted that these com- 
pleteness estimates are entirely based on point sources; for 
real galaxies with extended flux distributions, the complete- 
ness at a given magnitude will be lower. However, for the rel- 
atively bright galaxies considered in this analysis (K < 22.4, 
or 1 magnitude brighter than the worst-case 95% complete- 
ness limit), this should not be an issue. Also, this does not 
take into account the possibility of close pairs of objects be- 
ing blended, and thus mistakenly being counted as a single 
object by SExtractor; however, for this analysis such effects 
are likely to be small (further discussed in ® . 

4.6. Catalog format 



FIG. 6. — Comparison of our photometric redshifts to 119 spectroscopic 
redshifts in the UDS field (top panel) and fractional err or (bottom pane l ). Red 
filled circles denote those points flagged as QSOs in the Simpson et al] 420061) 
sample, and green open circles are objects with log(x") > 2.9 as derived by 
EAZY. 

The final ^-selected catalog employed in this analysis was 
generated from the 99022 objects (not lying near bad pixels) 
detected by SExtractor in the UDS mosaic. Some of these 
may be image artifacts and/or fall on bad or missing regions 
in the SXDS and SWIRE images as described in 503] Such 
objects are included in the catalog, but with flag(s) noting 
that their photometry in a given band is likely to be unreli- 
able. Photometry in the B through IRAC 4.5 /im bands, as 
well as limited morphological information, are included, and 
all fluxes are given with an AB zeropoint of 25. The columns 
of the catalog 11 are as follows: 

Columns (l)-(3). — Running ID, right ascension, and declina- 
tion (J2000) 

Columns (4)-(15). — BRi'z'JK "color" fluxes and errors (1 ."75 
diameter aperture, listed in the order fs, as, fa, &r, ■ ■■) 
Column (16). — Total K flux in the SExtractor AUTO aperture 
Columns (17)-(20). — IRAC 3.6/im and 4.5/im fluxes and er- 
rors (matched to the optical/IR "color" aperture) 
Columns (21)-(23). — ^f-band half-light radius, ellipticity, 
and position angle 

Columns (24)-(25).— Optical and SWIRE bad-pixel flags 
Column (26). — Flag for objects not detected in the stacked 
B+R+i'+z' image 

Column (27). — Internal flag generated by SExtractor 

5. PHOTOMETRIC REDSHIFTS AND REST-FRAME COLORS 

5.1. EAZY fitting 

Photometric redshifts we re determined with t he EAZY 
code 12 , described in detail in Brammer et al. (2008). In its de- 
fault configuration, EAZY uses x 2 minimization to fit linear 
combinations of six basis templates to broadband galaxy spec- 
tral energy distributions; a /T-band luminosity prior and esti- 
mates of systematic errors due to template mismatch are also 
taken into account. These default settings have been demon- 
strated to provide reliable photometric redshifts for other K- 

Available from 

|http : / /www . strw . leidenuniv . nl/galaxyevolution/UDsl 

12 Code and documentation are available at 
|http : //www, astro . yale ■ edu/eazy| 
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selected samples (see Brammer et al. 2 0081) . 

To test the reliability of the UDS z p hot values we com- 
pared them to galaxies in this field with known spectro- 
scopic redshifts. A query using the NASA Extragalactic 
Database (NED) finds 96 published redshifts, most of which 
are old passively-ev olving galaxies at 0.8 < z < 1.2 from 
lYamada et alj {2005). An additio nal 60 redshifts of radio 
sources are from the catalog of [Simpson et al. (20061). Af- 
ter cross-correlating these objects with objects detected in the 
UDS K band, 119 spectroscopic redshifts remained; EAZY 
was able to find z p hot solutions for 1 10 of these (and found no 
X 2 minimum for the other 9). The left-hand panel of Figure [6] 
shows a comparison of the photometric and spe ctroscopic red- 
shifts of these galaxies. Objects flagged by Simps on et al.l 
(2006) as "QSO" or "XQSO" (X-ray emitting QSO), which 
may not be well-fit by typical galaxy templates if the AGN 
light contributes significantly to the overall flux, are shown 
as filled red circles in the plot. Of the aforementioned nine 
objects for which EAZY failed to find z p hot solutions, four of 
them fell in this AGN category. 

While the agreement between the photometric and spectro- 
scopic redshifts is quite good (particularly at z ~ 1), roughly 
8% of the points are outliers with fractional errors of 0.2 or 
more (and about half of these in turn are serious outliers, 
with |Az| /(l +z sp ec) > 0.5). Determining a priori from pho- 
tometry alone which objects are likely to be outliers is dif- 
ficult - although most of the QSOs are indeed strong out- 
liers, weak (or even strong) AGN activity in galaxies is not 
always evident. However, the \ 2 value returned by EAZY 
appears to be a somewhat reliable indicator: all objects with 
log% 2 > 2.9 have (1 +z p hot) values that deviate more than 20% 
from their measured spectroscopic redshifts. When objects 
with these high \ 2 values are removed from the sample, the 
remaining normalized median absolute deviation (NMAD) in 
Az/(1 +z spec ) is 0.033 with a median offset of Az = -0.013. 

At z ~ 1.2 the photometric redshifts appear to exhibit a 
small systematic offset of Az/(1 +z) ~ -0.05 (though due to 
the small number of points here it is difficult to assess the 
magnitud e of the offset). A sim ilar offset at the same redshift 
is seen by Brammer et al. (2008) even with more photometric 
bands included (notably U, V, and H), and the offset persists 
when different template sets and input parameters are used. 
Additionally, relatively small-scale "spikes" and "voids" (on 
the scale of Az = 0.1-0.2) in the photometric redshift distri- 
bution are apparent-see Figure|7] These may in part represent 
real variations in the large-scale galaxy distribution, but some 
are also likely to be numerical artifacts intrinsic to the photo- 
metric redshifts. However, since these effects are relatively 
small, and since in our analysis we employ redshift intervals 
(Az > 0.5) larger than the observed "spikes" in the z p hot dis- 
tribution, such effects are not likely to significantly affect our 
results. 

The rest-frame near- and mid-infrared properties of typi- 
cal galaxies are not as well constrained as the optical SEDs, 
so the galaxy templates employed by EAZY tend to be un- 
certain in this regime. This uncertainty is taken into ac- 
count during the z p hot calculations; nevertheless, the inclu- 
sion of the IRAC data could still in principle introduce sys- 
tematic effects. To check this, we re-ran EAZY with only 
the BRi'z'JK photometry and compared the results to those 
with the IRAC fluxes included. While the resulting z p hot val- 
ues were similar, the scatter and systematic offset increased 
substantially, with NMAD [Az/(1 +z spe c)] = 0.044 and me- 




12 3 4 

Z phot 



FIG. 7. — Distribution of apparent K magnitude as a function of photometric 
redshift (top panel), and number of galaxies as a function of Zphot (bottom 
panel). 

dian Az = -0.017. Moreover, although the number of ma- 
jor outliers is comparable, such outliers are no longer easily 
rejected: even with the x 2 cut an d excluding objects where so- 
lutions were not found, only three of the AGNs are excluded, 
and overall the worst outliers are no longer confined to the fits 
with the highest \~ values. 

Including IRAC data in the calculation of photometric red- 
shifts thus reduces the z p hot scatter and allows outlying points 
(including many AGNs) to be rejected from the sample . This 
is perhaps not surprising: for example. IStern e t al. (2005|) find 
that AGNs typically have redder IRAC [3.6] -[4.5] colors 
than typical galaxies. Since the galaxy templates employed 
by EAZY do not take possible AGN components into ac- 
count, the 4.5 fim excess in galaxies with nuclear activity in 
turn cannot be matched well by any of these templates. In- 
deed, most of the objects that either had high \ 2 values or no 
solutions from EAZY appeared to exhibit redder [3.6] - [4.5] 
colors than those galaxies with good z p h t fits. 

5.2. Monte Carlo analysis 

During the process of computing z p h ot values, EAZY also 
provides estimates on the redshift uncertainty and a probabil- 
ity distribution p(z) for each object. The p(z) curves in par- 
ticular are useful for estimating the redshift distribution of a 
large sample of objects, since the true distribution is likely to 
be broader than a histogram of the best-fit z p hot values (e.g., 
in an extreme case with a sample of z p hot = 2.00 galaxies, the 
true distribution would be much broader than a delta func- 
tion). To assess the consistency of the p(z) distributions, we 
performed 120 Monte Carlo iterations varying the input pho- 
tometry assuming Gaussian errors. Though this is a relatively 
small number of iterations, it gives a rough estimate of the 
z p hot uncertainty for any given object, and an accurate compu- 
tation of the redshift distribution of a large sample of objects. 

For each object, the median and dispersion among the 120 
Monte Carlo runs were determined and compared to the es- 
timates from EAZY. The median values for the UDS sam- 
ple were consistent with the best-fit values derived by EAZY 
from the original (unperturbed) catalog, as were the lcr uncer- 
tainty estimates. Furthermore, for subsets of the catalog, the 
redshift distributions derived from the Monte Carlo analysis 
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FIG. 8. — Rest-frame color-color diagram for A^b < 22.4, Zphot < 2.5 
UDS galaxies (black points). Red points denote spectroscopically-confirmed 
Z ~ 1 old passively-evolving galaxies (OPEGs) from Yamada et al] 420051) 
with li ttle or no detected line emission (Wa(OII) < 5 A). Bruzual & Chariot 
12003) evolutionary tracks of passively-evolving (red line) and constantly 
star-forming (cyan line) stellar populations from 0.1 - lOGyr are overplot- 
ted, and the arrow shows the effect of 1 magnitude of dust extinction. The 
extended star-forming track and the "quiescent clump" (which overlaps with 
the OPEG sample) are clearly separated. 

closely reflected the EAZY-derived p(z) distributions. The 
uncertainties derived by EAZY thus appear to be internally 
consistent for the UDS data. 

5.3. Interpolating rest-frame photometry 

At z ~ 2.5, the U, V, and J rest-frame bands fall roughly 
into the observed J, K, and IRAC 4.5/im bands respectively. 
Up to this redshift, intrinsic UVJ fluxes can thus be inter- 
polated from the observed data. Filter response curves for 
the observed bands, taking into account atmospheric absorp- 
tion and detector quantum efficiency, were downloaded from 
each of the WFCAM, Suprime-Cam, and IRAC websites. For 
rest-frame filter definitions we used the standard filter defi- 
nitions witho ut atmo s pheric absorption or detector response 
included: the lBesselll (119901) U and V curves, and the Mauna 
Kea definition ( Tokunaga et al. 2002) for J. Res t-frame fluxes 
were t hen interpolated following the method of R udnick et alj 
(2003) 13 . Uncertainties on the rest-frame fluxes were derived 
using the perturbed input catalogs and output redshifts from 
the Monte Carlo analysis described in 35.21 Note that these 
color uncertainties only take into account photometric errors 
(and the resulting z p h t errors), and do not include the uncer- 
tainties intrinsic to the templates used forz p hot fitting. 

6. EVIDENCE FOR QUIESCENT GALAXIES TO Z ~ 2.5 

6.1. The rest-frame colors of quiescent and star-forming 
galaxies 

We use the galaxy catalog derived from the 
UDS/SXDS/SWIRE data, along with the photometric 
redshifts and interpolated rest-frame colors described in 
the previous section, to analyze the rest-frame color distri- 
bution of galaxies out to z = 2.5. Beyond this redshift the 

13 This was carried out using the Inter- 
Rest script (E. Taylor et al., in preparation) ; see 
|http : //www . strw . leidenuniv . nl/~ent/InterRest| 



rest-frame J band begins to "fall off" the reddest filter in 
our catalog (IRAC 4.5/im), and rest-frame colors at z > 2.5 
are therefore less reliable. Stars are selected (and removed 
from the sample) using the two criteria: J - K < and 
(J-K) < Q.2(i' - A")-0.16. This two-color cut was derived 
by inspection of the J-K vs. i' — K color-color diagram, 
wherein stars form a tight, well-defined track. Additionally, 
we exclude all objects that fall on bad pixels or missing data 
regions in any band, as well as those with bad photometric 
redshift solutions (logx 2 > 2.9 as derived in ^5j, and apply a 
magnitude limit of K < 22.4. When all these criteria are met, 
the subsample analyzed hereafter contains 30108 galaxies 
between < z < 2.5. By comparison with 119 spectroscopic 
redshifts in the field we find typical photometric redshift 
errors of Az/(l+z) ~ 0.033, though this is almost entirely 
measured with galaxies at z < 1 .2; at higher redshifts the z p h t 
uncertainty is likely to be substantially larger. 

Figure |8] shows the rest-frame U — V vs. V — J (hereafter 
UVJ) diagram for this subset of UDS galaxies. A striking bi- 
modality emerges: one diagonal track extends from blue to 
red V—J, while a localized clump that is r ed inU-V but blue 
in V — J lies above this track. Previously, lLabbe et"ail (120051) 
and IWuyts et"aTI (120071) found that actively star-forming and 
quiescent galaxies segregate themselves in this plane, with the 
star-forming galaxies forming a diagonal track and quiescent 
galaxies populating mostly the upper left-hand region; with 
the far larger number of sources in the UDS, it is now ev- 
ident that the galaxies form a truly bimodal distribution in 
this plane. Indeed, this interpretation is supported here by 
both data and models: the red points overlying the "quies- 
cent clump" in Figure [8] are spectr oscopically-confirm ed old 
passively-evolving galaxies from Yam ada et al.l (12005), while 
the star-forming and q uiescent loci are reasonabl y coincident 
with the corresponding B ruzual & CharloJ ((2003) stellar pop- 
ulation models; a detailed analysis of the star-formation prop- 
erties of galaxies in the UVJ plane will be presented in a forth- 
coming paper (Labbe et al., in preparation). 

Essentially, the UVJ diagram allows the degeneracy be- 
tween red star-forming and red quiescent galaxies to be bro- 
ken: while galaxies with blue U — V colors in general ex- 
hibit relatively unobscured star formation activity, red galax- 
ies could be either quiescent galaxies with evolved stellar 
populations or dust-obscured starbursts. But since dust-free 
quiescent galaxies are blue in V -J, they occupy a locus in 
the UVJ plane that is distinct from the star-forming galax- 
ies, allowing the two populations to be empirically separated. 
Clearly, such a separation using a single color (such as U—V) 
would be fraught with problems: at best, the quiescent sam- 
ple derived in this manner would be contaminated by red star- 
bursts, and if the number of such starbursts were sufficiently 
high the bimodality would no longer be visible (see also Fig- 
ure [Toll. Note that although the U—B color is better at dis- 
tinguishing between a narrow break (characteristic of an old 
stellar population) and broader dust reddening when spectro- 
scopic redshifts are available, the larger uncertainties on pho- 
tometric redshifts do not allow sufficiently accurate U—B col- 
ors to be estimated. 

Even when photometric redshifts are derived using alterna- 
tive template set s within EAZY, or an entirely differ ent z p h t 
code (HYPERZ; lBolzonella. Miralles. & Pelldll2000h . the ba- 
sic bimodal shape of Figure [8] persists, indicating that this 
method is robust to the specific numerical technique em- 
ployed. 
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FIG. 9. — Rest-frame U — V vs. V—J colors of UDS K < 22 A galaxies in five redshift bins. The greyscale represents the density of points in the central region 
of each plot, while outlying points areplotted individually. The solid lines show the adopted divisions between the star-forming and quiescent galaxy samples 
at each redshift, defined in Equation [4] Median uncertainties in the rest-frame colors of quiescent and star-forming galaxies (only taking random errors into 
account) are shown in the upper left and lower right of each frame, respectively. 




FIG. 10. — Left panel: Number of red galaxies as a function of the normal distance from the diagonal dividing lines shown in Figure|9] the normal distance is 
defined such that, in these histograms, star-forming galaxies are to the left (i.e. negative values) and quiescent galaxies are on the right. Bimodality along this 
diagonal line is clearly evident up to z = 2. Right panel: Color-magnitude diagrams in the same bins; red and black points denote quiescent and star-forming 
galaxies (selected on the basis of the cuts in Figure |9) respectively. The green line is a fit to the red points in the 0.0 < z < 0.5 bin, and the inset plots show 
galaxy counts as a function of the distance in U - V from this line (with quiescent galaxies plotted separately as a red histogram). Note that no U - V bimodality 
is evident at z > 1.5, even though it is clearly seen in the UVJ plane. 
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6.2. Color Evolution 

Changes in rest-frame colors with redshift reflect the intrin- 
sic evolution of stellar populations, since all spectra have in 
principle been transformed to the same reference frame (bar- 
ring possible Zphot systematic errors). UVJ diagrams at five 
different redshifts (in bins of width Az p hot = 0.5) are shown 
in Figure [9] The most immediately apparent feature is that 
the observed bimodal distribution of star-forming and quies- 
cent galaxies is clearly seen up to z ~ 2. Although the scatter 
increases substantially at higher redshifts (to the point where 
it likely washes out any intrinsic bimodality at z > 2), most 
likely due to a combination of larger photometric redshift un- 
certainties and weaker observed fluxes, the two distinct popu- 
lations are nonetheless visible. 

It is also notable that the shape of the star-forming sequence 
appears to change: at the lowest redshifts the sequence curves, 
but this may be due to the effects of using a small aperture on 
relatively nearby galaxies (i.e., the outermost, bluer parts of 
galaxies with large angular sizes falling outside the "color" 
aperture). At z > 0.5, where the angular size-redshift rela- 
tion begins to flatten, the dust sequence indeed becomes more 
linear. The increase in star formation activity at higher red- 
shifts is also apparent in two ways: the fractional number of 
dusty, star-forming galaxies (the upper-right portion of the 
dust sequence) increases, and the entire dust sequence appears 
to move to redder V—J colors. Note that at 2.0 < z p hot < 2.5 
a concentration at red V — J and U — V appears; this mainly 
comprises galaxies with a substantial contribution from the 
red Ay = 2.75 template included with EAZY. This single tem- 
plate substantially improves the z p h t solutions of very dusty 
galaxies, but can lead to the observed discrete clump in color- 
color space. 

With increasing redshift the dead clump appears to steadily 
move to bluer U — V (with the median U — V changing by 
~ 0.15 magnitudes between the z = 0-0.5 and z = 1.5-2 
bins), while the V—J color remains more or less unchanged. 
This is expected from passive evolution of the stellar popula- 
tions in the clump, but relatively small systematic errors in the 
photometric redshifts (on the order of Az/(1 +z) ~ 0.05) can 
produce similar offsets in the average U — V color of quies- 
cent galaxies at z ~ 1.5. The z p hot estimates show a systematic 
offset that is smaller than this, about (z p hot _ Zsrjec)/(l +z spe c) ~ 
-0.01 over the full redshift range (see Figure|6}. Nonetheless, 
given the lack of spectroscopic redshifts at z > 1.3, we cannot 
ascertain whether there are larger systematic offsets at higher 
redshift, and therefore cannot accurately measure the color 
evolution of galaxies in the quiescent clump or absolutely de- 
fine the boundaries of the box in which they live. 

6.3. Sample Division 

Since the color bimodality is visible up to z ~ 2, we em- 
ploy empirical criteria to divide this sample into quiescent 
and star-forming subsamples. In each redshift bin of Fig- 
ure [9] where the bimodality could be seen, an initial diagonal 
cut was first made between the two populations. Histograms 
of red galaxy counts relative to this line were then derived. 
The position of the diagonal cut was fine-tuned (keeping the 
same slope at all redshifts) to fall roughly between the two 
peaks (except in the lowest-redshift bin, where the shapes of 
the quiescent and star-forming galaxy loci made this impossi- 
ble); Figure [TO] (left panel) shows these histograms along the 
adjusted diagonal lines. 

The adopted diagonal selection criteria for quiescent galax- 
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FIG. 11. — Rest frame UVJ diagram in four redshift ranges with bins 
color coded by the median specific star-formation rate of the galaxies within 
each bin; outliers in the UVJ plane are plotted individually. The z < 0.5 bin 
is not included here because the adopted SSFR proxies are not reliable at 
low redshifts. Quiescent galaxies lie in approximately the same region at all 
redshifts. 

ies between < z < 2.0 bin are: 

(U-V) > 0.88 x (V-/) + 0.69 [0.0 < z < 0.5] 
(U-V)> 0.88 x(V-/) + 0.59 [0.5<z<1.0] (4) 
(U-V) > 0.88 x (V-/) + 0.49 [1.0 < z < 2.0] 

Additional criteria of U — V > 1.3 and V—J < 1.6 are applied 
to the quiescent galaxies at all redshifts to prevent contami- 
nation from unobscured and dusty star-forming galaxies, re- 
spectively. The samples of star-forming galaxies are then de- 
fined by everything falling outside this box (but within the 
color range plotted in Figure|9| such that the very small num- 
ber of extreme color outliers are not included in either sam- 
ple). The two distinct populations are no longer visible at 
z p hot > 2.0, but the same dividing line used at z = 1—2 is shown 
for reference. Note that the exact placement of this division 
does not significantly affect the results presented herein (see 
also §7 for a more detailed discussion). Typical random color 
uncertainties for the star-forming and quiescent samples at 
each redshift interval are shown in the lower right and upper 
left corners, respectively, of each panel in Figure [9] 

6.4. MIPS 24/im measurements 

Between z ~ 1-2.5, the 24/im flux of galaxies is strongly 
correlated with dust-obscured star formation due to the pres- 
ence of redshifted 6- 12/im p olycyclic aroma tic hydrocarbon 
(PAH) features in this band dYan et al.ll200H e.g.). The ob- 
served 24/im fluxes of the UDS galaxies can thus be used to 
empirically confirm whether the UVJ criteria indeed select 
star-forming and quiescent galaxies. Since the K band, which 
at these redshifts corresponds to rest-frame optical/NIR light, 
is a rough tracer of stellar mass, the lA^imlK flux ratio 
then provides an estimate of the dust-obscured specific star- 
formation rate (SSFR). Likewise, the UV-to-^T flux ratio pro- 
vides an approximate measure of the unobscured SSFR, and 
thus the combination of these three bands provides a reason- 
ably robust estimate of total SSFR. 

Using the FIREWORKS catalog of IWuyts et all (l2008h we 
derive proxies for the IR and UV SSFR based on UV/K and 
24-nm/K flux ratios, and use these proxies to obtain order- 
of-magnitude estimates of the SSFRs of UDS galaxies; the 
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FIG. 12.— Rest-frame colors of bright (K < 22.4), 2.0 < z < 2.5 galaxies 
(top panel). Red galaxies with 1.3 < (U -V) < 2.0 are selected by the boxes 
shown, and the median 24/im MIPS flux for the galaxies within each box is 
plotted in the lower panel. The 24pm flux rises at r ed V —J, confirming the 
trend suggested in the corresponding panel of Figure II II 



details of this derivation are given in the Appendix. Fig- 
ure QT] shows the UVJ diagram in four redshift bins between 
0.5 < z < 2.5 (the bin at z < 0.5 is not shown, since the SSFR 
proxies are not expected to be effective at such low redshifts). 
Since the 24/im data are too shallow for most galaxies in this 
sample to be detected individually, we combine the SSFRs 
of multiple galaxies as a function of UVJ color as follows: 
at each redshift galaxies are divided into square UVJ color 
bins. If > 10 objects fall into a given bin, a square is plot- 
ted with a color corresponding to the median SSFR of the 
underlying galaxies (in essence "binning" the galaxies' SSFR 
estimates); otherwise, individual points are plotted. A strong 
trend is clearly seen from z = 0.5-2: the "quiescent clump" 
on average exhibits low SSFR, and is distinctly offset from 
the adjacent red star-forming galaxies. 

In the 2.0 < z < 2.5 bin of Figure|9]the scatter in the UVJ di- 
agram is too large for any intrinsic bimodality to be seen, and 
the trend in SSFR also appears weaker (Figure [TTJ, perhaps 
due in part to the smaller number of points at high redshift and 
passive evolution. At z > 2, we employ another approach to 
verify the presence of galaxies with relatively low SFRs: the 
top panel of Figure [12] again shows the z = 2-2.5 UVJ plot, 
with red galaxies selected by the rectangular boxes shown. 
The median 24/im flux of the N galaxies within each box, 
measured in 8" diameter apertures, is then calculated. To esti- 
mate the uncertainties in the median fluxes, N apertures were 
placed randomly on the MIPS image and the median 24/im 
flux within these apertures calculated. This was repeated 500 
times for each value of N, and the clipped standard deviation 
of these 500 trials was taken as the uncertainty in the me- 
dian of the galaxy fluxes for every bin of N galaxies. Note 
that since the randomly-placed apertures may contain bright 
and/or confusing 24/im sources, this error estimate inherently 
takes both background and confusion noise into account. 

The median 24/im fluxes within each bin are shown in 
the lower panel of Figure Q~2j indeed, 24/im flux increases 
strongly with V -J color, indicating that V -J does indeed 
provide a good proxy for star formation rate at z = 2-2.5 
(with a factor of ~ 5 difference in 24/im flux between the 



FIG. 13. — Rest-frame V magnitude distributions of the quiescent (red 
histogram) and star-for ming (black histogram) galaxies in the four redshift 
ranges shown in Figure 1101 The vertical dashed line shows the value of 
M*, at the central re dshift of each bin, assuming the z = 0. 1 SDSS value of 
IBlanton et all 420031) and evolving as M v (z) = My(0)-z- At high redshifts, 
quiescent galaxies compose roughly half of the total bright galaxy population, 
and a larger fraction at lower redshifts. 

bluest and reddest galaxies in V -J). The trends in 24/mi 
flux and SSFR strongly suggest that the UVJ selection criteria 
indeed separate quiescent from star-forming galaxies (even 
given the limited number of photometric bands available in 
UDS) at z < 2, and may work reasonably well up to higher 
redshifts, z < 2.5. Since the color bimodality can no longer 
be seen at z > 2 it is possible that the quiescent galaxy sam- 
ple at these redshifts contains somewhat more contamination 
from red, star-forming galaxies; furthermore, the rest-frame 
J flux estimate relies on the relatively shallow 3.6/4.5 /mi data 
at z > 2, which may result in subtle incompleteness effects in 
the sample at these redshifts. Therefore, for the remainder of 
our analysis we focus only on galaxies at z < 2 where the bi- 
modality can obviously be seen. It is also important to note 
that the presence of AGN may produce enhanced 24/im emis- 
sion and redder V -J colors, and thus mimic the photometric 
effects of star formation. As described in $5] most objects 
with known strong AGN activity do not have valid z p h t solu- 
tions and have likely been removed from the sample; even so, 
some weak AGN activity may remain and contribute to the 
observed correlation. 

6.5. Luminosities 

Absolute rest-frame V magnitude histograms of the num- 
ber counts of quiescent and star-forming galaxies in each red- 
shift interval are shown in Figure Qj] T he value of Mt at the 
center of each redshift bin, assuming the Bla nton et alj d2003l) 
SDSS value at z = 0.1 and evolving as M£(z) = M£(0)-z, is 
also shown. At bright magnitudes (My < -22.5) between z ~ 
1-2, quiescent galaxies contribute approximately as much 
to the galaxy number density as those th at are active l y star- 
forming. This is similar to the result of iKriek et alJ (2006), 
who found star formation rates consistent with zero in 9/20 
bright galaxies at z = 2.0-2.7, though their spectroscopic de- 
termination of SFRs is likely to produce fundamentally differ- 
ent "quiescent" samples than our rest-frame color definition. 
At lower redshifts (0.5 < z < 1 .0), the luminous galaxy popu- 
lation is comprised of an larger fraction of quiescent galaxies, 
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while galaxies with fainter V luminosities (and hence which 
typically have lower stellar masses) are dominated by those 
undergoing star formation. 

7. CLUSTERING OF QUIESCENT AND STAR-FORMING GALAXIES 

To investigate the relation between halo mass and star for- 
mation activity at z > 1, we calculate the clustering of the 
star-forming and quiescent galaxies (as defined by their rest- 
frame UVJ colors) at these redshifts following the method of 
Ouad ri et all ([2007a, 2008); a brief summary of this method 
follows. The angular correlation funct i on of e ach sample was 
computed using the lLandv & Szalavl d 19931) estimator, tak- 
ing into account all "bad" data regions (described in detail 
in 34.31 ). Figure [14] shows these functions for the quiescent 
and star-forming samples, with the former clearly exhibiting 
stronger clustering (in angular space) than the latter. 

Because of the large number of objects in this catalog at 
z = 1 - 2, it is possible to further break each of the quiescent 
and star-forming samples into bright and faint subsamples, 
thereby also investigating the effect of luminosity on cluster- 
ing strength. Angular correlation functions were recomputed 
for star-forming and quiescent galaxies within two luminosity 
bins divided at My = -22.0. This is approximately the mini- 
mum luminosity at which the z = 1 - 2 sample is complete — 
essentially all galaxies with My < -22 at z = 2 have observed 
K fluxes brighter than the K < 22.4 cutoff, while at M v > -22 
some galaxies (most of which are low-mass and blue) fall be- 
low the K flux limit. However, this should not affect the relia- 
bility of the rq determination itself since the Monte Carlo z p h t 
simulations ( 35. 2b fully account for differences in the redshift 
distributions. For each subsample we fit power laws to the an- 
gular correlation functions over the range 60"- 300", and esti- 
mated ro using the Limber projection (see lOuadri et al.[2 007a. 
for details). The lower cutoff ensures that we are only measur- 
ing the large-scale correlation function, and are not unduly in- 
fluenced by the clustering of galaxies within individual halos. 
The upper cutoff was chosen to min imize the importan ce of 
the integral constraint correction (see O uadri et al.l2008l) . The 
slope of the spatial correlation function is fixed to 1 .6, which 
is consistent with each of the subsamples studied here, and 
the uncertainties are estimated using bootstrap resampling. 

The resulting clustering lengths of the quiescent and star- 
forming samples in both luminosity bins are shown in Fig- 
ure[l5] The clustering of bright quiescent galaxies is stronger 
(at 1.4cr significance) than that of the bright star-forming 
sample, with r = 9.5+\ A 2 h" 1 Mpc and 7.4 ± 0.9 h" 1 Mpc re- 
spectively, and the difference in clustering length becomes 
far more pronounced in the fainter bin. 14 Interestingly, the 
clustering length of quiescent galaxies appears to be indepen- 
dent of luminosity; star-forming galaxies, on the other hand, 
exhibit a marginal increase in clustering strength with lumi- 
nosity. The stronger clustering of bright galaxies in the com- 
bined (star-forming plus quiescent) sample thus appears to 
be driven by both the larger quiescent fraction among bright 
galaxies as well as the stronger clustering of star-forming 
galaxies. 

The exact criteria used to divide the "quiescent" and "star- 
forming" samples (i.e. the lines shown in Figure|9]l may have 
an effect on this result. While we chose the diagonal dividing 

14 Since the objects which are "missed" in the faint luminosity bin are pre- 
dominantly low-mass blue galaxies, the true clustering of the star-forming 
sample in this bin is likely to be even lower than what we find; i.e. the dif- 
ference between ro for the faint quiescent and star-forming galaxies would 
increase with a deeper K flux limit. 
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FIG. 14. — Angular correlation functions for the 1.0 < z p hot < 2.0 quies- 
cent (red filled circles and errorbars) and star-forming (blue open circles ) 
samples. The quiescent galaxy sample shows significantly stronger angular 
clustering. 

line to lie centered between the star-forming track and qui- 
escent "clump," this is not necessarily the best criterion. For 
example, moving the line "upward" (toward redder U — V and 
bluer V -J) would result in a less complete sample of quies- 
cent galaxies that is also less contaminated by star-forming 
galaxies scattering across the line; similarly, moving the line 
down toward the star-forming track increases completeness at 
the expense of contamination. To check whether this affects 
the clustering measurement, we moved the diagonal divider 
by 0.1 magnitudes in both directions and re-calculated the ro 
values. The perturbed ro values were all within 1 a of the old 
values, with one exception: when the dividing line was moved 
upwards (i.e. defining a less complete quiescent galaxy sam- 
ple with less contamination), the clustering strength of bright 
(My < -22) quiescent galaxies increased by ~ 1.5er. There- 
fore, the result that quiescent galaxies exhibit a larger cluster- 
ing length appears robust to the color criteria used to select 
the quiescent sample. 

The clustering length we have found for quiescent galaxies 
at z = 1-2 is fully consistent with the value o f rp = 10.6 ± 
1-6/T 1 Mpc for Zphot = 2-3 DRGs in this catalog (Oua driet alj 
2008). This is perhaps not surprising, as roughly half of the 
Zphot > 2 DRGs in this sample fall within the "quiescent" UVJ 
selection region. 

8. DISCUSSION 

The presence of a bimodal population of galaxies in the 
UVJ plane to at least z ~ 2 is consistent with p revious spec- 
troscopic studies, such as iRriek et all (l2008bl) . who find a 
nascent red sequence at z ~ 2.3. The purely photometric anal- 
ysis presented herein allows much larger samples of quiescent 
galaxies to be found, albeit with lower confidence for any in- 
dividual object. At first glance, howeve r, this appears to be in 
conflict with the Cirasuol o et al.l d2007l) claim that the galaxy 
bimodality disappears at z ~ 1.25- 1.5. In fact, the apparent 
discrepancy bet ween these results is probably due to observa- 
tional effects: Cirasuolo et al. (2007) consider only bimodal- 
ity in the rest-frame color-magnitude (U-B vs. Mb) diagram, 
i.e. the "red sequence" method. Indeed, as the right panel of 
Figure [10] shows, such a plot does not show a bimodal popu- 
lation of galaxies at z p hot > 1-5, even though the bimodality is 
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FIG. 1 5 . — Deprojected real-space clustering length in two absolute magni- 
tude bins (My < —22 and My > —22) for the 1 .0 < Zphoi < 2.0 quiescent (red 
filled circles and errorbars), star-forming (blue open circles), and combined 
(black crosses) galaxy samples. For each point, the median My of the corre- 
sponding sample is plotted. The lower panel shows the fraction of quiescent 
galaxies in each bin. Quiescent galaxies exhibit stronger clustering, and thus 
occupy more massive dark matter halos, than star-forming galaxies in both 
magnitude bins. 

still present in the UVJ diagram. Similarly, the absence of a 
bimodal distribution in the UVJ diagram at z > 2 (Figure [TOl 
left panel) does not necessarily imply the disappearance of 
bimodality in the underlying galaxy population, but rather re- 
flects the larger photometric and z p hot uncertainties at these 
redshifts (as can also be seen from the large UVJ error bars at 
2 < Zphot < 2.5; Figure |9). 

Using a single rest-frame color (determined via photomet- 
ric redshifts) to select quiescent and star-forming galaxies can 
obscure an underlying bimodality for two reasons: (1) the 
population of dusty, star-forming galaxies increases with red- 
shift, and these have significant overlap their optical colors 
with any red, quiescent galaxies that may be present; and (2) 
the greater scatter in the rest-frame U - V or U - B at higher 
redshifts can effectively wash out the intrinsic red sequence. 
This increased scatter is due not only to the larger uncer- 
tainties on measured photometry of fainter galaxies at large 
redshift, but also th e larger photometric redshift uncertainties 
dKriek et al.ll2008ab . As previously mentioned, U — B by it- 
self is only effective at distinguishing quiescent from dusty 
galaxies when accurate spectroscopic redshifts are available; 
with photometric redshifts the rest-frame color uncertainty 
is too large to see the quiescent red sequence. The use of 
V—J in addition to U — V helps mitigate both of these prob- 
lems: V — J separates red, quiescent galaxies from red, dust- 
obscured starbursts (Figure [SJ, and a two-dimensional color 
bimodality is more resilient to increased scatter than a sin- 
gle color. Thus, when only photometric redshifts are avail- 
able, the two-color UVJ technique employed herein appears 
to provide a "cleaner" selection of high-redshift, quiescent 
galaxies than standard red-sequence methods. 

The differences between high-redshift quiescent and star- 
forming galaxy properties can provide clues to the physical 
mechanisms behind star formation (or the lack thereof). In the 
moderately high-redshift (z = 1-2) sample considered here, 
the most striking contrasts between the two populations are in 
their clustering strengths and luminosities. Quiescent galax- 



ies cluster more strongly than those undergoing active star 
formation (implying, from standard CDM theory, that they 
reside preferentially in higher-mass halos) and typically ex- 
hibit brighter absolute V magnitudes, comprising about half 
of the bright galaxy population at all redshifts. However, the 
clustering strength of quiescent galaxies does not appear to 
be strongly dependent on luminosity: as shown in Figure [151 
the bright and faint quiescent samples exhibit roughly equal 
clustering lengths, which in turn are consistently higher than 
the corresponding values for star-forming galaxies. This may 
be indicative of a characteristic halo mass above which star 
formation is inhibited. 

Our results may be taken to suggest that the star formation- 
density relation was already in place at z ~ 1 .5, with quiescent 
galaxies more strongly clustered and hence situated, on aver- 
age, in denser regions. On first glance this would seem to 
contradict recent claims of a "reversal" of the star formation- 
densi t y relation that occu rs at z ~ 1 — 1.5 (e.g. Elb az et ail 
l2007t ICooper et all 120081) . These authors present evidence 
that the average SFR per galaxy in dense regions is higher 
than in less dense regions, in contrast to the well-known re- 
lationship in the local universe. Whether this constitutes a 
"reversal" is in part a matter of terminology; it may be true 
that the typical star formation rate increases with en viron- 
mental density (as found by, e.g., ICooper et al.l 120081) . but 
this does not preclude the possibility that the relatively rare 
quiescent galaxies are preferentially found in the densest re- 
gions (as suggested by our results). We also note that the star 
formation-density relation may be qualitatively different de- 
pending on whether one is considering a bsolute star formatio n 
rates or specific star formation rates (see ICooper et al.l 2008). 

The observed segregation of quiescent and star-forming 
galaxies both in clustering strength and luminosity suggests 
strong links between halo mass, stellar mass, and the ces- 
sation of star formation. Different theoretical bases for the 
quenching of star formation in galaxies hosted by massive 
halos have been proposed, some focusing on radio-mode 
AGN feedback proces ses which prevent further gas cooling 
and star formation (e.g. lCroton et al.ll2006l:lBower et al.ll2006i 
and references therein), while others suggest that hot accre- 
tion and shock processes in high-mass halos may be suffi- 
cient to produce quiescent galaxies without invoking AGN 
jKeres et all [20051; iBirnboim. Dekel. & Neisteinl 120071) . Al- 
though these models may not represent the definitive expla- 
nation for the observational results presented here, the corre- 
spondence between theory and observation is nonetheless in- 
triguing. Further observations of large quiescent galaxy sam- 
ples (e.g., comparing the radio emission of high-z quiescent 
and star-forming galaxies) are likely to shed further light on 
the underlying quenching mechanism. 

9. CAVEATS 

This analysis employs photometric redshifts, which of 
course are much less well-constrained than those determined 
using spectroscopy. Furthermore, the number of spectro- 
scopic redshifts for comparison is quite small, with 119 avail- 
able spec-zs in this field. Essentially all of these known red- 
shifts are from two samples, neither of which can be consid- 
ered unbiased: a set of z ~ 1 old passively-evolving galax- 
ies, and galaxies selected from their radio emission (thus all 
likely exhibiting some degree of AGN activity). If anything, 
one might expect photometric redshift determinations to be 
worse in this latter sample since the templates used to de- 
rive the phot-zs do not include AGN activity; however, Fig- 
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ure [6] shows that the fits are actually quite good from about 
z = 0- 1.5 (and the worst outliers, including two-thirds those 
with signatures of AGN, are excluded through either a x 2 cut 
or the lack of a z p h t solution). The fit to the "old" z ~ 1 
galaxy sample was even better, with fractional redshift errors 
on the order of 2%. At z > 1 .5 it is likely that these fractional 
errors are at least ~ 0.07, the scatter found in the EAZY- 
derived photometric redshifts from the CDFS-GOODS cata- 
log (which included several m ore photometric ban ds than the 
UDS+SXDS+SWIRE dataset; iBrammer et alj|2008b . 

It is nonetheless possible that some additional systematic 
errors remain in the photometric redshifts. Such errors would 
manifest themselves as shifts in the rest-frame colors, the 
magnitude of which primarily depend on the galaxy redshift 
and template shape. However, the results presented herein 
are not dependent on the exact rest-frame colors themselves, 
but rather on the relative separation between the quiescent 
and star-forming populations. Since the two populations are 
clearly visible from z = 0-2 (Figure |9), and the positions of 
galaxies in the UVJ plane are strongly correlated with 24/im 
flux (FigurefTTTi. we conclude that we can distinguish popula- 
tions of galaxies based on their relative star-formation rates. 
It is nonetheless important to note that the absolute values of 
the interpolated rest-frame colors may still be subject to sys- 
tematic, possibly redshift-dependent offsets, and these off- 
sets may be different depending on the SED (e.g. for quies- 
cent and star-forming galaxies). Further differences are likely 
to come about with differing photometric redshift codes and 
techniques for determining rest-frame colors. This under- 
scores the importance of defining such empirical color cuts 
based on the actual data and analysis techniques employed; 
relying on criteria defined on other datasets can lead to inac- 
curate results. 

Another possible source of error is the ability of SExtrac- 
tor to deblend close pairs of objects. Where two extended 
sources overlap, or if the separation of two point sources is 
smaller than a few times the K-band PSF FWHM (~0"7), 
oftentimes the two objects will be counted as a single object. 
In this case, the galaxy clustering is probably underestimated 
at the smallest scales (< 1 .5" or so). Even if they are success- 
fully identified as separate objects, the fixed color apertures 
used to measure fluxes are likely to be contaminated, and thus 
may cause errors in photometric redshifts and rest-frame col- 
ors. However, the fraction of objects in such close pairs is 
negligible compared to the total catalog. Furthermore, these 
clustering results are based on larger-scale correlations — i.e., 
we do not calculate the clustering below 9 ~ 2". Thus, this 
incompleteness at the smallest scales is unlikely to affect our 
results. 

10. CONCLUSIONS 

From our matched multiband UDS+SXDS+SWIRE K- 
selected catalog we have derived reliable photometric red- 
shifts and rest-frame colors. For a relatively bright (K < 22.4) 
subsample of this catalog, we find that: 

1 . Galaxies show strong bimodal behavior in rest-frame 
U — V vs. V—J color-color space, with one population 



of "dead" (non-star forming) galaxies, and a sequence 
of dusty, actively star-forming galaxies. This behav- 
ior can be seen out to z ~ 2, so the observed present- 
day galaxy bimodality was present at least to this red- 
shift; however, at z > 1.5 the bimodality is not seen in 
a single-color-magnitude diagram. 
2. From z = to 2.5, the rest-frame V—J color (at red U— 
V colors) shows a strong correlation with our estimates 
of specific star-formation rate, indicating that the V -J 
color is a reliable tracer of star formation activity even 
at z > 2 where the rest-frame color uncertainties are too 
large for the bimodality to be seen. 



3. Quiescent and star-forming galaxies at z > 1 contribute 
roughly equally to the overall galaxy population at the 
brightest V luminosities; the less luminous population 
is dominated by actively star-forming objects. 



4. The clustering strength of quiescent galaxies from z = 
1-2 appears to be independent of galaxy luminosity, 
and is consistently stronger than the clustering of the 
actively star-forming galaxies; thus, "dead" galaxies 
appear to occupy more massive halos than those which 
are in the process of forming stars. This suggests a link 
between halo mass and the early cessation of star for- 
mation activity. 

Forthcoming data releases of the UKIDSS UDS, including 
deep //-band data and ultimately reaching a depth of K = 24.9 
(AB), will allow these results to be tested with even fainter, 
higher-redshift galaxies. More importantly, the 290-hour 
ultra-deep Spitzer legacy survey of the UDS field (SpUDS; PI: 
J. Dunlop) currently being undertaken will provide at least an 
order-of-magnitude increase in exposure time compared to 
the SWIRE data used in this paper, giving much-improved 
constraints on photometric redshifts, AGN activity, and rest- 
frame colors. 
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APPENDIX 

ESTIMATING SPECIFIC STAR-FORMATION RATES 

Accurately determining specific star formation rates (SSFR; i.e. SFR as a fraction of galaxy mass) typically requires deep UV 
and 24/im data in order to measure unobscured and obscured star formation, respectively. Unfortunately, the UDS data presented 



16 



Williams et al. 




FIG. 16. — Relations used to derive specific star -formation rates fr om the observed IR (fi4um/fK', 'op panel) and UV (fit/ fx', foottom panel) flux ratios. Data 
in this figure are from the FIREWORKS catalog o flWuvts etal] 420081) ; the IR relation shown is fit to galaxies between 0.5 < z < 2.5 and the UV at 1 .0 < z < 2.5 
(a similar UV relation was derived for the z = 0.5 - 1 galaxies, not shown here). The IR SSFR shows a good correlation with fu^ m / '/k over the full range, but 
the UV relation tends to underestimate large SSFRs. 

here currently lack several key ingredients; i.e. there are no U -band data, the available SWIRE 24/im data are shallow, and the 
necessary detailed models for determining galaxy masses and the 24/im-L/R conversions are beyond the scope of this work. 
However, to test the validity of the UVJ colors as a tracer of SSFR, such precision is not necessary: simple order-of-magnitude 
estimates will suffice. 

Over the redshift range 0.5 < z < 2.5, we therefore res ort to a proxy based entirely on observed photometry and calibrated with 
the Chandra Deep Field-South FIREWORKS catalog of IWuvts et all d2008l) . Reliable SSFRs were derived from FIREWORKS 
based on the rest-frame 2800A luminosity and observed 24/im fluxes (Wuyts et al., in preparation). The sum of these two contri- 
butions should then approximate the total amount of star formation in "typical" galaxies (i.e. those without extreme obscuration, 
like SMGs). Since K flux is itself a tracer of stellar mass, it follows that the /uv //k A ux ratio (for some observed band that falls 
in or near the rest-frame ultraviolet) can provide a rough estimate of the unobscured SSFR, and /24 M „,/ 'fx should do the same for 
the "dusty" SSFR. 

From the FIREWORKS data we find that a single linear relation sufficiently describes the correlation between the IR SSFR 
and fiA^mj '/k over 0.5 < z < 2.5. There is also a strong correlation between UV SSFR and fa/ fx', however, the exact slope 
and normalization of the relation appear to change abruptly at z ~ 1 ; for the UV relation we thus perform separate linear fits at 
0.5 < z < 1 and 1 < z < 2.5. This UV proxy is somewhat flatter (i.e. insensitive to fit/ fx) than the 24/im-SFR relation, and 
substantially underestimates the largest UV SSFRs by a factor of ~ 2, but it appears to accurately reproduce relatively low SSFRs 
(la2 x 10 9 yr _1 ); thus, this UV proxy is effectively used as an additive term to correct the fluxes of galaxies which are faint in 
24/im. The best-fit relations, shown in Figure[l6] are: 

SSFRir = -2.2x 10- n + 1.85 x 10- 10 (/24///r) 

SSFR uv (z< l) = -7.4x 10- 10 +3.8x \Q-\f R /f K ) (Al) 
SSFR uv (z> l) = -5.8 x 10" 10 +5.8 x 10" 9 (/«//^) 

Note that the UV relations formally yield substantially negative SSFR values even when fu/fx > 0, since SSFR UV is not quite 
a linear function of /r/Zk', thus, we require that the UV contribution to the total SSFR be positive (i.e. for any object where 
SSFRuv < 0, we set SSFRuv = 0). Since a large number of galaxies have /24 M „, < due to background noise, and the SSFRir 
relation appears very close to linear, we do allow individual objects to have SSFRir < in the stacking analysis. 
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